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ABSTRACT 

We study possible anomalous CP-violating Higgs couplings to /i/i and ti fully 
model-independent way through top-quark pair productions at muon colliders. As- 
suming additional non-standard neutral Higgs bosons, whose couplings with top- 
quark and muon are expressed in the most general covariant form, we carry out 
analyses of effects which they are expected to produce via CP- violating asymmetries 
and also the optimal-observable (00) procedure under longitudinal and transverse 
muon polarizations. We find the measurement of the asymmetry for longitudinal 
beam polarization could be useful to catch some signal of CP violation, and an 00 
analysis might also be useful if we could reduce the number of unknown parame- 
ters with a help of other experiments and if the size of the parameters is at least 
0(1) ~ 0(10). 
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1. Introduction 

It is widely known that the standard model of the electroweak interaction (SM) has 
been so far quite successful in describing various phenomena below the electroweak 
scale with high precision. Its top-quark and Higgs-boson sectors are however still 
not fully-tested part of the model. If there exists any new physics beyond the SM 
within our reach, its effects will be likely to appear in those sectors. Therefore it 
is worth to look for experiments which allow for a comprehensive investigation of 
top-quark and Higgs-boson properties. 

Anomalous top-quark interactions could be tested, for instance, at the ee col- 
liders in the International Linear Collider (ILC) project p,j. However it is not easy 
to study Higgs sector thereby. Muon colliders were proposed as an ideal machine 
to explore Higgs properties [2j. From a purely theoretical point of view, muon 
colliders are quite similar to ee colliders, but the fact that a muon is much heavier 
than an electron could provide with a non-negligible difference in phenomenological 
studies of Higgs sector. 

Indeed many authors have studied how to analyze Higgs-top interactions at 
muon colliders. Most of them focused on the resonance region, i.e., direct Higgs 
productions, and/or —>■ (Higgs) ti in the framework of some specific models 
with multi Higgs doublets, like MSSM, and pointed out that a muon collider will 
be a useful tool to identify CP properties of Higgs scalars |2]-[1]. 

As a complementary work to them, we study in this article possible anomalous 
Higgs interactions with and ti in a fully model-independent way through 
ti processes. Our main purpose is to clarify to what extent we would be able to draw 
a general conclusion on those interactions without assuming any particular models 
at muon colliders in off-resonance region. In other words, we aim to study the 
possibility and limit of muon colliders for model-independent analyses of possible 
new physics in the top-quark and Higgs-boson sectors. 

After describing our calculational framework in section 2, we compute CP- 
violating asymmetries for both longitudinal and transverse beam polarizations 
in section 3, where based on the results we also discuss a detectability of the 
anomalous-coupling parameters. In section 4, we study whether the optimal- 



observable procedure is effective when we try to determine the anomalous param- 
eters separately. Finally, a summary is given in section 5. 

2. Framework 

As mentioned in Introduction, we perform a model-independent analysis of pos- 
sible non-standard Higgs interactions with top-quark/muon for longitudinal and 
transverse beam polarizations. Let us summarize our framework first which is the 
basis of our calculations. Throughout this paper, we express the standard-model 
Higgs as h and the non-standard neutral Higgs as H. 
Effective amplitude 

The invariant amplitude of /i/i — > (7, Z, h, H) — > ti corresponding to Figure [1] 
is given as follows: 



M{^lJi tt) = M^i + Mz + Mh + Mh, 




where A^7,z are the standard 7 and Z exchange terms 
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Figure 1: Feynman diagrams for fifi —>■ (7, Z, h, H) —>■ tt 



D-y^zis) being the propagators multiplied by the coupling constants 

D,{s) ^ -^e^i , (6) 

''^^'^"lecos^^^.-Ml' 
e, g, 9w being the elementary charge, the SU{2) gauge coupling, the Weinberg 
angle respectively, and s = {p/j, + Pp)^, 
Mhis the standard Higgs-boson exchange term 

Mh = Dh{s) u{pMPt) ■ v{p-^Mp^), (8) 

while A4h is the non-standard Higgs exchange contribution, for which we assume 
the most general covariant form: 

Mh = Dh{s) u{Pt){at + ha^)v{pi) 

X + btilb)u{p^), (9) 



with 



As = 10 

mf — s — imdi 



(i = H), Fi and v being the total decay width and the vacuum expectation 
value of the SM Higgs field. We treat a^^^ and bt^j^ as complex numbers to take into 
account the possibility that they are form factors. 

Readers may claim that we assume only one additional Higgs-boson in spite of 
our statement that we perform a fully model-independent analysis. In fact, our 
frame can incorporate any number of Higgs exchange terms. It will be clear by 
re-expressing such an amplitude as 

N 

Al[Non-SM Higgs] = J^DhAs) uiPt){ai + 6*75)^ (pf) ■ v{p^^{a^^ + 1f^-f5)u{p^) 
1=1 

^J24^ii^Hi{s) u{Pt)v{Pt) ■ v{Pfi)u{p^) 

i 

i 
i 

+ EKKDHA^)^(PthAPt)-v{Pi,Hu{p^)- (11) 



Thus, all the contributions can be packed into our parameters as follows: 

ata^ = Y.a^^a'^DHXs)/DHAs). (12) 

i 

atb, = J2a%DHXs)/DHAs), (13) 

i 

ha^ = Y.blai,DHXs)/DHAs), (14) 

i 

b,b, = J2b%DHAs)/DnAs). (15) 

i 

Beam polarization 

The beam polarization, P, along one axis (polarization axis) whose direction is 
defined by a unit vector s is given by 

p=P±l^P^^ (16) 

where p±s is the number density of the particle in each beam whose spin component 
on this axis is ±s. We can take into account this polarization by multiplying the 
spin vector s° in the projection operator u{p)u{p) and v{p)v{p) by P. That is, we 
are to use (0, Ps) as the spin vector in its rest frame. 

In the following, we choose the direction of p^ as the z axis and express the 
azimuthal angle of s as (p. Then the degree of the longitudinal polarization is 



given by Pl = Psz, that of the transverse polarization by Py = yP^ — P|, and 
consequently 

(0,Ps) = (0, Pt cos 0,Pt sin 0,Pl). (17) 

The /i and p, spin vectors in the fifi CM frame are obtained from f[T7|) via appropriate 
Lorentz transformations as 

= [PliP, Pt cos (/), Pt sin Pli) , (18) 
s° = {PliI3, Pt cos 0, Pt sin 0, -Pli) , (19) 

where 



/? = y^l - Aml/s, 7=1/^1-/?' (20) 
and the momenta of pi and /2 in this frame are 

p" = ^v^(l,0,0,/3), p° = ^v^(l, 0,0,-/3). (21) 



3. CP-violating asymmetries 

It is straightforward to calculate the cross section a{fifi ti) starting from ampli- 
tude ([1]) as 

^ -aififi tt) = ^ l'^^\ \M{fifi tt)\'^ . (22) 



dcos9 

We perform this via FORM [5] , but the analytical result is a bit too long to give here 
explicitly. Therefore we show in the following our results numerically. Throughout 
our analysis in this article, we take \Pl\ = 1 or \Pt\ = 1. It may seem to be an 
extreme and unrealistic assumption, but we chose the polarization this way because 
our aim here is to know "to what extent" we could know about the anomalous 
interaction (Q, i.e., we would like to study the possibility and limit of muon colliders 
for model-independent analyses of new physics beyond the standard model. 
Numerical results 

We study two CP-violating asymmetries and A^, the former of which is the 
one for longitudinal beam polarization 

a{++) + a{ ) 

and the latter is the one for transverse polarization 

a(x = Jr/2) + ff(x = -T/2)' 

where cr(±±) express the cross sections for = = ±1, while a{x = ±7r/2) are 
the ones for P^ = Pt = I with x = 4' — 4> = ±vr/2. We here chose |x| to be tt/2 
since it maximizes the CP- violation effects (see, e.g., [4]). 

Concerning the decay widths of h and H, Ph^Ui they are of course different 
quantities, but we use the same formula for Ph as Ph within the standard model 
[6] (see the later discussions). The other SM parameters are taken as follows: 

sin^ Ow = 0.23, Mz = 91.187 GeV, v = 246 GeV, 

mt = 174 GeV, = 105.658 MeV, ruh = 150 GeV. 

In Figures [2] and [3] are presented A^ as functions of y/s and mn, while in Figures 
m and [5] are given At in the same way for Reaj_^ = Imat,^ = Rebt,^ = Imfej^^ = 0.2 




Figure 3: mn dependence oi Al 
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Figure 4: ^/s dependence of At 
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as an example to sketch a rough feature of these quantities. We find that the 
absolute value of could be sizable, but that of is very small. 

It may seem strange that there appears such big difference between \Al\ and 
\At\ though they are both CP- violating asymmetries computed with the same 
anomalous parameters. The reason is in their denominators. In the case of Al, 

not only the numerator cr(++) — cr( ) but also the denominator receives little 

contribution from 7/Z exchange terms, while they can contribute to cr(x = ±7r/2) 
without being suppressed except in the difference a{x = ^/2) — cr(x = — 7r/2). 
Indeed if we focus on the numerators alone, there is only small difference between 
Al and At- For ^/s = 550 GeV and rriH = 500 GeV with the same anomalous 
couplings as in the figures, e.g., 

Al ■■ - or( ) = 1.6 X 10~^ fb, (25) 

At : a{x = 7r/2) - a{x = -7^/2) = 4.8 x 10"^ fb. (26) 

It is obvious that the peaks in A^^t are all due to the H propagator, but readers 
might wonder why At changes its sign while Al not in the vicinity of s = Mfj. 
Therefore, it would also be helpful to give a brief explanation here about those 
different behaviors of them. As mentioned above, a{x = ±7r/2) receive sizable 
contributions from 7/Z exchange terms. This means the interference between the 
7/Z-exchange and if-exchange terms, which is proportional the H propagator, is 
important in At, and its sign changes thereby depending on whether s > Mfj or 
s < Mfj. On the other hand, this is no longer the case for Al since the 'j/Z- 

exchange terms are suppressed in both a{++) and a{ ). Therefore the sign of 

Al is determined by the difference of cr(±±), where the size of the amplitude A4h 
itself is much more crucial. Here we chose = 6^ as an illustration, which makes 

cr(++) larger than cr( ) and leads to positive Al as will be understood from 

eq. (j34l) on //(++) and a similar calculation for A4h{ )• This also tells us that 

different parameters, e.g., = —b^ could make Al negative. 
Detectability of the asymmetry 

Let us study the detectability of Al, that is, the expected statistical precision 
in its measurement, which tells us how precisely we would be able to determine 



Al. For instance if we take y/s = uih =500 GeV with Rea^^^ = Refe^^^ = Imot^^ = 
Im&f^ = 0.2, Al becomes 0.73, while the cross sections are (t(++) = 5.0 x 10~^ fb 

and cr( ) = 7.8 x 10"'^ fb, leading to ~ 29e events for an integrated luminosity 

L = 500 fb^^, where we expressed the detection efficiency of tt productions as e. 
They are combined to give the following statistical uncertainty: 

SAl = ^(1 - Al)/N = 0.68/VeL = 0.13/v^. (27) 

Consequently, the expected statistical significance Nsd is 

NsD = \Al\/SAl = 5.7v^. (28) 

That is, we can confirm l^^^l 7^ at 5.7y^ level. For example, Nsd = 4.0 for 
e = 0.5. Here, assuming L = 500 fb~^ may be a bit too optimistic, but we used 
this value considering that we aim to find the possibility and limit of the muon 
colliders as mentioned in the beginning of this section. It is easy to transform our 
numerical results for any other L. 

We have given an example of Nsd for = rriH =500 GeV, but it is not 
general, so let us show the results for some other ^/s in Table [Tj changing also the 
parameters as Reaj,^ = Re6f_^ = Imaf_^ = Imbt^^ = 0.1, 0.2, 0.3. There we find 
that we would be able to observe some signal of CP violation as long as we are not 
too far from the H pole. 
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Table 1: Nsd as a function of y/s (with e = 1 for simplicity) for Reat^^ = Rebt^^ = 
Imaf_^ = Im&t^^ = 0.1 (a), 0.2 (b), and 0.3 (c) 

We used the SM formula for Fh as an appropriate approximation (/^ =67.5 
GeV for rriH =500 GeV [6J), since we did not introduce any new light particles 



that can appear in the final state of h and H decays. Strictly speaking, however, 
a new mode like H ^ hh might be possible for = 150 GeV and Mh = 500 
GeV. Instead of re-computing Fh including such new modes, which demands us to 
assume a concrete form of those couplings, we give Nsd for = 80 and 100 GeV 
with ReOf^^ = Refe^^^ = Ima^^^ = Imbt^^ = 0.2 in Tabled which tells us that our 
conclusion would not be affected so much, especially in the off-resonance region. 

^/s (GeV) Fh = Fhimn) Fh = 80 GeV Fh = 100 GeV 
450 1.5 1.4 1.3 
500 5.7 4.4 3.2 
550 2.4 2.3 2.1 
600 LO U) 1.0 

Table 2: Nsd as a function of ^/s for F^ = 80 and 100 GeV, where F^ = Fh{mii) 
means that Fh was computed with the SM formula {Fh{mH) = 67.5 GeV). 

Parameter dependence of the asymmetry 

Measuring would be quite interesting, but what receives is of course one 
single combination of contributions from all the anomalous parameters. We are 
performing a model-independent analysis of possible new-physics effects, but once 
we get actual experimental data, our results are going to be applied for a realistic 
model building, li Ai does not depend on some parameters so much, it will be hard 
to test any models in which those parameters play a significant role. Therefore it 
must be important to see how A^ depends on each parameter. 

Let us study how Nsd changes when we vary one parameter from 0.0 to 0.3. 
The results are given in Tables [3] and HI where Nsd are presented for one of the 
parameters = 0.0, 0.1, 0.2 and 0.3 with the others being fixed to be 0.2. There 





Re at 


Re 


Re6i 


Re 


Im at 


Im 


Im ht 


Imfe^ 


0.0 


4.1 


4.4 


4.7 


2.2 


4.1 


4.4 


4.7 


2.2 


0.1 


4.7 


5.3 


4.9 


3.8 


4.7 


5.3 


4.9 


3.8 


0.2 


5.7 


5.7 


5.7 


5.7 


5.7 


5.7 


5.7 


5.7 


0.3 


7.0 


5.6 


7.1 


7.7 


7.0 


5.6 


7.1 


7.7 



Table 3: Nsd as a function of each parameter for ^/s = 500 GeV with the rest 
being fixed to be 0.2 





Re at 


Re 


Rebt 


Re b^ 


Im at 


Im 


Im bt 


Imfe^ 


0.0 


1.1 


2.0 


2.0 


0.4 


2.3 


1.5 


1.9 


1.8 


0.1 


1.7 


2.3 


2.1 


1.3 


2.3 


1.9 


2.0 


2.1 


0.2 


2.4 


2.4 


2.4 


2.4 


2.4 


2.4 


2.4 


2.4 


0.3 


3.2 


2.4 


2.9 


3.6 


2.6 


2.8 


2.9 


2.5 



Table 4: Nsd as a function of each parameter for ^/s = 550 GeV with the rest 
being fixed to be 0.2 



we observe that Nsd receives a contribution from every parameter though there 
are some differences among them, which indicates that any model will be testable 
through measuring Al. 

4. Optimal-observable analysis 

The optimal-observable technique [7] is a useful tool for estimating expected sta- 
tistical uncertainties in various coupling measurements. Suppose we have a cross 
section 

^(=r(0)) = ^cj,(0), (29) 

where are known functions of the final-state variables and q's are model- 

dependent coefficients. The goal is to determine the q's. This can be done by using 
appropriate weighting functions Wi{(j)) such that / Wi{(j))E{(j))d(j) = Ci. In general 
different choices for Wi{(j)) are possible, but there is a unique choice for which the 
resultant statistical error is minimized. Such functions are given by 



w. 



:(0) = ^X,,/,(0)/r(0), (30) 



where Xij is the inverse matrix of Aiij which is defined as 

M^, = j^^dct>. (31) 

When we use these weighting functions, the statistical uncertainty of Cj is obtained 
as 



5ci = ^X,,aT/N, (32) 
where ax = J{da/d(j))d(j) and is the total number of events. 



We study whether we could get more information of the anomalous parameters 
via this procedure. Here we focus on the longitudinal beam polarization, since 
we found that it is practically impossible to catch any new-physics signal for the 
transverse beam polarization even when we could fully use the total cross sections. 

In oder to apply this technique to our analysis, we need to express the an- 
gular distribution of the produced top quark in terms of the anomalous-coupling 
parameters like eq. (l29|) . We have altogether eight independent parameters since 
we assumed the all couplings a^^^ and 6^,/^ to be complex. Although our aim is to 
perform an analysis as model-independently as possible, it will be too complicated 
to treat them all equally. Therefore we here assume that the size of the imaginary 
part of each parameter is much smaller than that of its real part. Since the imagi- 
nary part of parameters (form factors) is often produced through higher order loop 
corrections in an underlying theory, this assumption is not unreasonable. We also 
drop the terms quartic in the anomalous parameters. 

With this reduced parameter set and assumption, the top-quark angular distri- 
bution should be represented as 

-r^a++{fifi tt) 
dcosO 

= fsuiO) + CaafaaiO) + CabfabiO) + Cbafba{0) + Cbbfbb{0), (33) 

where fsuiO) expresses the SM contribution, 

Caa = (Re at) (Re a^), Cab = {Re at) {Re bf,), 
Cba = (Re6t)(Rea^), Cbb = {Rebt){Reb^), 

and fij{0) {i,j = a,b) are all independent of each other. If we reverse the signs of 

Cab and Cba, we get da {^fi — > tt). 

Practically, however, fia{0) and fib{0) become equivalent in the limit of 0. 
We can see this as follows: In calculating da+j^, the muon-spinor part becomes 

1 + 75 

+ bi^l5)u+{p^) - v+{pf,){a^ + b^-f5)^—u{p^) 

= (a^ + bt.)v+{Pf,)u+{p^) (34) 
in the limit. That is, a^^ and 6^^ contribute almost equally to da++, which leads to 

faa{9) ^ fab{9), fba{e) ^ fbb{e). 



Since we keep finite, it is in principle possible to perform an analysis treating 
all fij{0) as independent functions, but it is clear that we end up having very poor 
precision thereby ^ . Therefore we neglect their differences from the beginning and 
start from 

ti) ^ M9) + Caf2{e) + c,fs{e), (35) 

dcost^ 

where = Caa + Cab, Cb = Cba + Cbfe, fi{0) = fsuiO), f2{d) = faa{d) ~ fba{d) and 

W) = fbaie) ^ ue). 

Using these functions we obtain the following results as the matrix fIST]) for 
^ = 550 GeV and mu = 500 GeV: 

Mu = 7.75 ■ 10"^ Mi2 = 5.80 ■ 10"^ M13 = -2.38 ■ 10"^ 

M22 = 4.35 ■ 10~\ M23 = -1.79 • 10-2, M33 = 7.37 ■ 10-^ (36) 

where we used fi{9) for S{(f)) in eq. fl3T]) . We then compute the (2,2) and (3,3) 
elements of the inverse matrix of M: 

X22 = 3.68 ■ 10^ X33 = 5.38 ■ 10^ . 

This means the expected statistical uncertainty in Ca^b measurements are 

5ca = 1.92- 10^ /^/L, 5cb = 2.32- 10^ /Vl. (37) 

This tells us that we need L = 3.7-10^ fb^^ for achieving 6ca = 1 and L = 5.4-10^ 
fb^^ for 6cb = 1, which are both far beyond our reach! 

We then assume one of the parameters is determined in some other experiments 
in order to look for realistic solutions. First, if Cq was unknown (i.e., if Cb was 
measured elsewhere), the corresponding precision becomes 6ca = 44.5/-\/L, i.e., 
6ca = 1.99 for L = 500 fb^^. We give also results for some other ^/s in Table [5l 
Conversely, if Cb is undetermined (i.e., only Ca is known), we have 5cb = 539/\/L, 
i.e., 6cb = 24.1 for L = 500 fb^^. Some other results are in Table O 

Therefore, if the size of Cq is 0(1), there is some hope to catch new-physics 
signal thereby. On the other hand, \cb\ is required to be at least 0(10). Note that 
it is never unrealistic to assume \ca^b\ to be 0(1) ~ 0(10) as is known in various 
models with two (or multi) Higgs-doublets (see, e.g., |9] and the references therein). 
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Table 5: Expected precision of Ca determination for niH = 500 GeV 



(GeV) 


Scb 


400 


168//L 


450 


227//L 


480 


329//L 


520 


447/ /L 


580 


878//L 


600 


im/VZ 



Table 6: Expected precision of Cb determination for niH = 500 GeV 



What we could know via Al measurements is only on CP violation, while Ca,b 
are both combinations of CP-conserving and CP-violating parameters. Therefore 
those two approaches could work complementarily to each other. 

5. Summary 

We have carried out a model-independent analysis of possible non-standard Higgs 
interactions with tt and through top-quark pair productions at future muon 
colliders. As was pointed out in Refs.[2]-[1], the muon colliders are quite useful for 
studying the Higgs sector around the resonance. Considering those preceding stud- 
ies, our main purpose here was to see if we could also draw any useful information 
in the off-resonance region without depending on any specific models. 

Starting from the most general covariant amplitude, we computed two CP- 
violating asymmetries for longitudinal and transverse beam polarizations in order 
to see if we could get any signal of new-physics which breaks CP symmetry, and also 
studied whether we could determine the non-standard-coupling parameters sepa- 
rately through the optimal-observable (00) procedure as a more detailed analysis. 



We found that the longitudinal CP- violating asymmetry Al would be sizable, 
while the transverse asymmetry At is too small to be a meaningful observable. We 
then estimated the detectability of ^4^ and showed that we would be able to observe 
some signal of CP violation as long as we are not too far from the H pole. We also 
studied there in some detail how Al depends on each parameter, and found that 
we have no parameter that contribute little, although there are some differences 
among the parameters. 

On the other hand, more detailed analyses via the 00 procedure seem challeng- 
ing. However, if we could reduce the number of unknown parameters with a help 
of other experiments, and if the size of the parameters is at least 0(1) ~ 0(10), we 
might be able to get some meaningful information thereby. Readers may claim that 
the use of the asymmetry Al is enough when we have only one unknown parameter, 
but this is not necessarily true. What we could draw from Al is information on 
pure CP violation, while we could also know something about CP-conserving part 
through an 00 analysis. 

In our approach, we need the total cross section of /ip, tt and the angular 
distribution of the final top quark, for which we only have to reconstruct the top- 
quark jet axis. If we further try to study, e.g., the final lepton distributions in the 
top-quark decays, we could get additional information on possible anomalous tbW 
coupling, however in that case we would suffer from another suppression factor, 
i.e., the branching ratio of the top-quark semileptonic decay. Therefore it will be 
more advantageous to use the top-quark distribution as a whole when performing 
off-resonance analyses at muon colliders. 

ACKNOWLEDGMENTS 

We are grateful to Federico von der Pahlen for giving us a valuable comment 
about the integrated luminosity we used in the text and also pointing out some 
typos in the manuscript. This work is supported in part by the Grant-in- Aid 
for Scientific Research No. 13135219 and No. 16540258 from the Japan Society for 
the Promotion of Science, and the Grant-in- Aid for Young Scientists No. 17740157 
from the Ministry of Education, Culture, Sports, Science and Technology of Japan. 



The algebraic calculations using FORM were carried out on the computer system 
at Yukawa Institute for Theoretical Physics (YITP), Kyoto University. 



REFERENCES 



[1] See ILC web site http : //www . linearcollider . org/ cms/J 



[2] V.D. Barger, M.S. Berger, J.F. Gunion and T. Han, Phys. Rev. Lett. 75 (1995), 



1462 (hep-ph/9504330); 



V.D. Barger, M.S. Berger, J.F. Gunion and T. Han, Phys. Kept. 286 (1997), 
1 ( Ihep-phyMoMTs] . 

[3] B. Grzadkowski and J.F. Gunion, Phys. Lett. B350 (1995), 218 



(|hep-ph/9501339j); 



D. Atwood and A. Soni, Phys. Rev. D52 (1995), 6271 ( |hep-ph/9505233| ); 
S.Y. Choi and J.S. Lee, Phys. Rev. D61 (2000), 111702 ( |hi^-ph /99d93l5|; 

E. Asakawa, A. Sugamoto and I. Watanabe, Eur. Phys. J. C17 (2000), 279 ; 
E. Asakawa, S.Y. Choi and J.S. Lee, Phys. Rev. D63 (2001), 015012 



(hep-ph/0005118). 



See also: A. Pilaftsis, Phys. Rev. Lett. 77 (1996), 4996 ( lhep-ph/9 603328|. 



[4] B. Grzadkowski, J.F. Gunion and J. Pliszka, Nucl. Phys. B583 (2000), 49 



( |hep-ph/0003091| ). 



[5] J. A.M. Vermaseren, Symbolic Manipulation with FORM" , version 2, Tutorial 
and Reference Manual, CAN, Amsterdam 1991, ISBN 90-74116-01-9. 

[6] A. Djouadi, J. Kalinowski and M. Spira, Comput. Phys. Commun. 108 (1998), 



56 (hep-ph/9704448) 



[7] D. Atwood and A. Soni, Phys. Rev. D45 (1992), 2405; 

M. Davier, L. Duflot, F. Le Diberder and A. Rouge, Phys. Lett. B306 (1993), 
411; 

M. Diehl and O. Nachtmann, Z. Phys. C62 (1994), 397; 



J.F. Gunion, B. Grzadkowski and X.G. He, Phys. Rev. Lett. 77 (1996), 5172 



(hep-ph/9605326). 



B. Grzadkowski, Z. Hioki, K. Ohkuma and J. Wudka, Nucl. Phys. 



B689 (2004), 108 (hep-ph/0310159); Phys. Lett. B593 (2004), 189 



(hep-ph/0403174); JHEP 0511 (2005), 029 (hep-ph/0508183) 



J.F. Gunion, H.E. Haber, G. Kane and S. Dawson, "T/ie Higgs Hunter's 
Guide" , Perseus Publishing, Cambridge, Massachusetts. 



